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Abstract. - We study the kinetics of domain growth of fluid mixtures quenched from a 
disordered to a lamellar phase. At low viscosities, in two dimensions, when hydrodynamic 
modes become important, dynamical scaling is verified in the form C{k,t) ~ L" f[{k — kM)L] 
where C is the structure factor with maximum at ItM and L is a typical length of the system with 
logarithmic growth at late times. The presence of extended defects can explain the behavior 
of L. Three-dimensional simulations confirm that diffuse grain boundaries inhibit complete 
ordering of lamellae. Applied shear flow alleviates frustration and gives power-law growth at 
all times. 



Introduction. - The kinetics of the growth of ordered phases after a quench from a disor- 
dered state continues to provide interesting physical questions. While the process is reasonably 
well understood for the case of binary mixtures, where dynamical scaling occurs and domains 
grow with power-law behavior [1], basic questions remain to be clarified for more complex 
systems [2]. In this Letter we consider the case of fluid mixtures where, due to competing 
interactions, the system would arrange itself in stripes. Examples are di-block copolymer 
melts, where chains of type A and B, covalently bonded end-to-end in pairs, segregate at 
low temperatures in regions separated by a stack of lamellae [3] , or ternary mixtures where 
surfactant form interfaces between oil and water [4]. Other systems with lamellar order are 
smectic liquid crystals [5], dipolar [6] and supercooled liquids [7], chemically reactive binary 
mixtures [8] . Lamellar patterns are also observed in Raleigh-Benard cells above the convective 
threshold [9]. 
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The ordering of lamellar systems is characterized by the presence of frustration on large 
scales. This affects the late time evolution which, as discussed below, can be very slow or 
also frozen. Relations with the dynamics of structural glasses have been also considered in 
recent literature [10, 11]. One can expect that in real systems the effects of the velocity field, 
inducing motion around local or extended defects, can be relevant. Our purpose here is to 
analyze the process of lamellar ordering taking into account the effects of hydrodynamics. 
We consider a model based on a Ginzburg-Landau free-energy, where dynamics is described 
by convection-diffusion and Navier-Stokes equations [12]. Our results for two-dimensional 
systems show that, when hydrodynamical modes are effective (at sufficiently low viscosities), 
dynamical scaling holds with a complex dependence of characteristic lengths on time, analo- 
gous to that of systems with quenched disorder [13] not present here. An intermediate faster 
growth is followed by a slower logarithmic evolution at later times. The slowing-down of the 
dynamics is attributed to the formation of extended defects between domains of lamellae with 
different orientation, almost perpendicular each other. Our simulations show that the same 
phenomenon occurs in three-dimensional systems. These results are not expected to depend 
on the specific model considered and could be relevant for a broad class of systems with 
lamellar order. We also verified that applied shear flows alleviate frustration giving power-law 
growth at all times. 

The model. - We consider the free-energy model [14, 15] 

TW} = J d'xC-^^ + + f I P +|(VV)n (1) 

where ip is the order parameter field representing the local concentration difference between 
the two components of the mixture. The parameters b and c have to be positive in order 
to ensure stability. For a < and k > the two homogeneous phases with = ±y''—a/b 
coexist. A negative k favors the presence of interfaces and a transition into a lamellar phase 
can occur. In single mode approximation, assuming a profile like Asinkox for the direction 
transversal to the lamellae, one finds the transition (|a| = 6) at a « — 1.11k^/c where fco = 
-y/— k/2c, = 4(1 -I- K^/4cb)/3. The expression (1) can be also written for negative k as 
J^{ip} = J d'^x{^(p'^ + jip'^ + f [(V^ -I- fco)'P)^]}' T — a — cfcg, and in this form is generally used 
to describe di-block copolymers in the weak segregation limit [15]. 
The dynamical equations are [2] 

—^+v-Vva, = + lyV^o. , 2 

ot g OX/3 

f + v.(.„,.rv»|^ , (3, 

where Va are the components of the velocity field, g is the total density of the mixture and 
the incompressibility condition V • t> = is considered, v is the kinematic viscosity and 
r is a mobility coefficient. The pressure tensor is the sum of the usual hydrodynamical 
part and of a tensor P^lf"^ depending on ip, with a functional form obtainable from the 
free energy containing off-diagonal terms, whose expression can be found in [12]. It can 
be shown that 5J- /S^p = implies daP^^^ = 0. The laplacian in the r.h.s. of Eq. (O 
ensures the conservation of the order parameter. We will simulate these equations by using a 
finite difference scheme for the convection-diffusion equation and a Lattice Boltzmann Method 
(LBM) for the Navier-Stokes equation. Advantages of this method [16] with respect to other 
LBM schemes for fluid mixtures are that spurious terms appearing in Eq.(3) [17] are avoided 
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Fig. 1 Fig. 2 

Fig. 1 - Configurations at different times of a 256 x 256 portion of the lattice. 

Fig. 2 - Time behavior of characteristic size L and Cm = C{kM,t) on log- log (inset) and log- linear 
scales. The straight line in the inset has slope 0.3. 



and the numerical cfRcicncy is increased. In _D = 3 we used a parallel version of the LBM code 
which fuses the streaming and the collision processes to reduce bandwidth requirements [18]. 

We can summarize now what is known on the ordering properties of lamellar fluids. 
Previous two-dimensional studies of the full model (1-3) [12], without quantitative analy- 
sis due to the small size of the lattice considered (128^), showed the relevance of hydro- 
dynamics for obtaining well ordered domains on the scale of the system simulated. More 
results exist for variants of the model (1) with long-range interactions, without hydrody- 
namics [19-21] or neglecting the inertial terms of the l.h.s. of Eq.(2) [22]. Eq.(3) without 
advection and with non-conserved order parameter, corresponding to the Swift-Hohenberg 
model for Raleigh-Benard convection [23 27], has been also largely studied. Regimes with 
dynamical scaling have been found with the order parameter structure factor behaving as 
C{k,t) =< |<pfc(t)p >~ i^/[(fc — fco)t^] with different values for the exponent z, ipk{t) being 
the Fourier transform of (p(x,t) [22,24-30]. Slower evolution with frozen states and grain 
boundary pinning has been also observed [10, 26] for very deep quenches. In £> = 3 the few 
existing simulations have not considered dynamical scaling [31,32]. 

Ordering dynamics. - We run our two-dimensional simulations on a 1024^ lattice starting 
from disordered configurations. We fixed \a\ = b so that the minimum of the potential part 
of the free-energy is at <^ = ±1. We checked on small systems (128^) for different k in the 
lamellar phase that the expected equilibrium state was reached. The value of c was fixed in 
such a way that the period is about 10 times the lattice spacing and |a| = 2 x 10~^,k = 
—3 X \a\,c = 3.8 X |a|,r = 25 for the cases presented. Details on the LBM part of the code 
can be found in [33]. 

After the quench, lamellae start to form evolving until the equilibrium wavelength is locally 
reached. This part of the ordering process does not depend on the value of viscosity. Later, 
however, the system continues to order only at sufficiently low viscosities {v ^ 0.1). Otherwise, 
local defects dominate [12] and, at times of order t ~ 3000, the system results frozen in tangled 
configurations. An example of evolution at low viscosity {ly = 8.33 x 10~^) is shown in Fig.l. 
Different kinds of defects can be observed. The annihilation of a couple of dislocations is 
put in evidence. We have measured the relative velocity of the dislocations and we found it 
constant. At the last time of the figure the system can be observed to be ordered on larger 
scales. 

A quantitative description of the ordering process comes from the analysis of the structure 

factor C{k,t). After spherical average we plotted C{k,t) at different times. In the early 

regime, with t ~ 3000, C{k, t) develops a maximum at a momentum kM which decreases with 
time until the equilibrium value ko is reached. Then the peak of C{k,t), remaining at fco, 
continues to grow while the width decreases indicating an increase of order in the system. A 
characteristic length L for this process can be extracted from the structure factor in the usual 
way by measuring the full width 6k of C at half maximum and defining L{t) = 2ir/dk. In the 
inset of Fig. 2 we plotted L(t) on log-log scale, calculated by averaging over 5 different runs. 
For more than one decade a behavior consistent with the power-law L ^ can be observed 
with z = 0.30 ± 0.02. The growth becomes slower at later times and we find i ~ Int as the 
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Fig. 3 Fig. 4 

Fig. 3 - Configuration of a 128'^ system at t — 1000000. The following parameters were used: 
a = -b = -0.026, K = -0.005, c = 0.0025, F = 0.25, and v = 0.1. 

Fig. 4 - The middle horizontal sections of a 512^ systems at times t = 600000 and t = 700000. 
Parameters of this simulations are the same of Fig. 3. Observe how the vertical stripes in the central 
bottom part of the figures are distorted by the pressure of horizontal domains. 



straight line in the main frame of Fig. 2 suggests. We checked that the slowing down cannot 
be attributed to finite size effects (at the latest times considered, L{t) is less than 1/10 of the 
size of the system) . 

The behavior of L can be related to the role played in late time dynamics by extended 
defects (grain boundaries) between domains of differently oriented lamellae. A decrease in 
the growing rate of L was observed at times t « 25000 when the average value of the angle 
between the normals of neighboring domains becomes close to 90°. Isolated grain boundaries 
can be shown to be stable defects in lamellar systems described by Eq.(^) [10]. A quantita- 
tive description of ordering in lamellar systems can be obtained by considering the evolution 
equation for a single defect. The position Xgb of the grain boundary evolves according to 
Xgb = AC^{t) — B cos{2koXgb + 4>) where C is the curvature of lamellae parallel to the boundary 
and A,B,(f> are model-dependent parameters [10]. If one assumes that dynamical scaling is 
verified, a typical length can be defined as the inverse of curvature, and the previous equa- 
tion implies a growth like t^^^ until a critical value for curvature is reached. This analysis, 
elaborated for one single defect, suggests that in the ordering of a large system, in absence 
of other processes and without hydrodynamics, defects will become pinned at a certain time 
with the system frozen in a configuration with many grain boundaries. In our case the velocity 
field helps the system in continuing the ordering process also at very late times but slower. 
We measured L ~ Ini which is the growth expected for systems with metastable states. An 
extension of the analytical arguments of [10,34] in presence of hydrodynamics is not available. 
However, it is interesting to observe that our simulations suggest an intermediate regime with 
a power law consistent with that of [10] followed by a logarithmic asymptotic growth, both 
results explainable as due to grain boundary dynamics. 

We also considered three dimensional systems. Simulations ranging from 64'^ to 512^ 
lattices were performed: The largest grid asked for about 20 days of sustained computations, 
using a 8-CPU vector machine NEC SX6 and more than 32 GB of memory. While complete 
order was reached on 64"^ lattices, larger systems confirm the role of grain boundaries in 
lamellar ordering. An example is given in Fig. 3 where a single grain boundary is shown: After 
t ^ 10^ no clear 2D ordering was observed. On the 512^ lattice we studied the behavior of 
L{t), defined as in _D = 2. After the formation of lamellae, the intermediate and late time 
regimes were found shorter and it was not possible to deduce quantitative behaviors. However, 
the role of grain boundaries can be appreciated looking for example at the configurations of 
Fig. 4. Here the central horizontal sections of the system at t = 600000 and t = 700000 are 
shown. Many grain boundaries with three main orientations can be seen. After t ~ 600000, 
in correspondence of a sudden decrease observed in the growing rate of L{t), the further 
evolution of the system proceeds through the distortion of the grain boundaries existing at 
that time. This can be seen for example in the central bottom part of the configuration at 
t — 700000. Different domains, trying to increase their size, push themselves reciprocally, 
helped by the velocity field. This gives an unexpected late-time increase of curvature, also 
observed in D — 2. The role of hydrodynamics is important in this process: The system 
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Fig. 5 Fig. 6 

Fig. 5 - Rescaled structure factors at different times. Symbols A, o, *, •, *, o refer to times 4600, 
14100, 40800, 61100, 78100, 123100, respectively. In the inset Cm is plotted on log-log scale as a 
function of L; the slope of the line is 1.25. 

Fig. 6 - Evolution of Lx and Ly for 7 = 10"''. The straight line has slope 0.6. At the latest times of 
the figure the system is almost completely ordered, as shown in the 200 x 200 portion of the inset. 



can reach a higher degree of order only by developing a strong velocity field which bends the 
existing flat grain boundaries. This will eventually bring later to a configuration with less 
defects. This picture is confirined by the increase of kinetic energy observed at f ~ 600000. 

Dynamical scaling. - In order to analyze dynamical scaling, together with L{t), we also 
considered the behavior of the peak Cm = C{kM, t) of the structure fac;tor. In two-dimensional 
simulations, from Fig. 2, one sees that Cm grows similarly to L{t) but with a different rate. 
This suggests to consider scaling behavior in the general form C{k,t) ~ L°'f[{k — kM)L] [35]. 
This differs from the usual expression considered in lamellar ordering (C ^ f[{k — kM)^^]) 
for the introduction of a and the use of L even if not given by a power law. As illustrated 
in the inset of Fig. 5, we evaluated a = 1.25. Then, by using this value, we found, as shown 
in Fig. 5, that rescaled structure factors overlap at all times after t ^ 3000 on a single curve, 
confirming our scaling assumption and suggesting an analogy with kinetic behavior in systems 
with quenched disorder. In Random Field Ising Models (RFIM), for example, diffusive growth 
changes over to logarithmic behavior but dynamical scaling holds at all times with the same 
scaling function as in standard Ising model [13]. Here a similar scenario occurs but without 
quenched disorder. This is interesting also in relation with the existence of an equilibrium 
glass transition in systems with lamellar order [11]. 

Concerning the exponent a, its value is related to the compactness of domains and is 
a = £> in asymptotic growth of binary mixtures, a = in early time regime when interfaces 
are forming [35]. Our results with a > 1 show that the argument for a = 1 based on the 
existence of one scaling dimension for 2D lamellar systems [24] does not hold in presence of 
hydrodynamics. We found dynamical scaling also in 3D simulations but for a shorter time 
interval. 

Shear effects. - We also studied lamellar ordering in presence of shear flow. Shear was 
applied with bounce-back conditions on rigid walls moving with velocity iz;^,^;; at the top and 
bottom of the simulation volume [33]. Ordering is favored since interfaces want to orientate 
with the flow [36]. Breaking and recombination of domains, induced by shear, makes the 
elimination of topological defects faster [37]. As a result, logarithmic growth disappears and 
ordering occurs also asymptotically with power-law behavior even in cases, at high viscosity, 
where without shear freezing would have been observed. 

An example of this behavior {v = 7.5) is shown in Fig. 6. Growth was measured by the 
quantities La = J dkC{k,t)/ J dk\ka\C{k,t) {a = x,y). Lx grows with exponent Zx = 0.6 
while Ly relaxes to the equilibrium value. Similar results were found for other values of u; a 
detailed exposition of our results with shear will appear elsewhere. 

Conclusions. - To conclude, we have studied lamellar ordering in fluid mixtures with 
competitive interactions. At low viscosities, when hydrodynamic modes are effective, dynami- 
cal scaling holds in the form C{k, t) ^ L" f[(k — kM)L] where L is a characteristic length with 
complex behavior and logarithmic growth at late times. This scaling, similar to that found in 
systems with quenched disorder, is shown here for the first time in 2D systems with lamellar 
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order. 3D simulations also show the relevance of grain boundaries defects for this dynamics. 
Shear flow removes extended defects giving power-law behavior at all times. 
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